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1.  INTRODUCTION 


An  important  element  in  our  National  Defense  is  the  development  of 
satellite  bassd  surveillance  systems  to  monitor  the  appearance  and  progress 
of  any  potential  threat.  The  Infrared  (IR)  spectral  region  for  detection 
systems  shows  good  promise  for  being  able  to  detect  various  targets  at 
different  altitudes.  In  support  of  continued  development  of  IR  detection 
systems,  it  is  necessary  to  fully  specify  and  understand  the  upwelling  IR 
radiance  from  the  earth-atmosphere  scene  within  a  given  field-of-view  (FOV) . 
Detecting  the  radiation  from  a  target  requires  that  not  only  the  target 
intensity  but  also  the  spectral  and  spatial  structure  of  the  background 
be  fully  understood.  Onboard  processing  capabilities  must  be  able  to  dis¬ 
tinguish  variations  in  the  upwelling  radiation  due  to  various  target  signatures 
from  variations  arising  from  background  spatial  and  temporal  features.  Thus, 
any  detection  system  must  be  designed  to  know  and  interpret  the  variances  in 
the  levels  of  IR  radiation  from  a  FOV  scene;  these  variations  include  the 
temporal,  spatial,  and  spectral  structure  of  the  radiation. 

Important  Inputs  to  the  design  of  successful  detection  systems  are  the 
magnitude  and  variations  in  the  radiance  from  various  sources.  This  knowledge 
will  aid  in  developing  data  processing  algorithms  to  distinguish  radiation 
signatures  of  clouds  and  other  surface  features  from  those  of  the  target. 
Spectral,  spatial,  and  temporal  data  combined  with  phenomenological  analyses 
will  provide  the  information  required  for  systems  design.  The  data  and 
analysis  will  address  such  questions  as  fluctuation  intensities,  signal  to 
noise,  solar  scattering  effects,  FOV  stability,  and  effects  of  atmospheric 
turbulence.  In  many  cases,  these  data  are  best  summarized  in  validated  cal- 
culatlonal  models  which  quantitatively  describe  the  important  physical  phenomena 
and  are  convenient  for  use  by  the  systems  designer. 


Clouds  will  be  a  significant  source  of  variations  in  the  IR  radiation 
for  any  system  viewing  the  earth’s  surface,  as  they  cover  approximately  .30 
percent  of  the  earth  surface  at  any  given  time.  Clouds  affect  a  satellite 
system's  ability  to  detect  a  target  from  space  in  several  different  ways: 

e  Obscuration:  Clouds  occult  the  sensor's  view  of  the  target. 

•  Background  Signal  Variations:  Clouds  change  the  background  radiance 
level  in  a  FOV  because  they  scatter  solar  radiation  and  have  a 
blackbody  temperature  different  from  that  of  the  earth  surface. 

e  Signal  Degradation:  Optically  thin  clouds  modulate  the  target  and 
background  signals  by  reducing  their  intensity.  Additionally, 
high  altitude  cirrus  will  also  scatter  solar  radiation  into  the 
line  of  sight. 

The  Balloon  Airborne  Mosaic  Measurements  (BAMM)  Program  is  tasked  to 
make  extensive  measurements  of  the  earth  surface  as  would  be  seen  from  a 
high  altitude,  down-looking  system.  Since  the  BAMM  platform  is  at  an 
altitude  of  30  km,  these  data  include  low  and  high  altitude  clouds,  as  well 
as  earth  surface  features  such  as  land-water  interfaces.  In  support  of 
this  program  Aerodyne  Research  has  developed  a  general  calculational  model 
for  the  upwelling  earth  radiance,  based  on  the  AFGL  LOVTRAN  code  for 
atmospheric  transmittance  and  radiation.  The  calculational  model  for 
scattering  and  emission  of  radiation  by  a  cloud  uses  an  approximate  solu¬ 
tion  to  the  radiative  transfer  equation,  which  incorporates  multiple  scat¬ 
tering  effects.  The  purpose  of  the  model  is  to  enable  quick  calculations 
of  the  upwelling  atmospheric  radiation  for  arbitrary  sun-viewer-cloud- 
background  viewing  scenarios.  The  primary  motivation  for  the  modeling 
approach  is  to  develop  a  calculationally  flexible  model  which  captures 
the  salient  features  of  the  IR  background  scene.  Therefore,  the  model  is 
based  on  the  LOWTRAN  atmospheric  transmission  code  and  a  simplified  diffu¬ 
sion  scattering  model  for  describing  multiple  scattering  effects  within  a 
cloud.  In  this  report  the  model  is  compared  to  radiance  data  for  clouds 
viewed  from  the  BAMM  and  other  measurement  platforms. 
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Section  2  of  the  report  discusses  modeling  of  the  physical  processes  for 
scattering,  absorption,  and  emission  of  radiation  within  a  cloud.  The 
important  assumptions  which  are  made  for  describing  the  radiation  transfer 
are  identified.  This  section  also  shows  the  application  of  a  simple 
two-stream  radiation  transfer  model  for  the  radiation  field  and  includes 
a  detailed  description  of  the  formalism  for  calculating  the  radiation 
from  a  cloud.  In  Section  3,  we  compare  the  model  predictions  to  data  taken 
by  the  Arthur  D.  Little  Company  and  by  the  Air  Force  Geophysics  Laboratory 
Flying  Laboratory  (KC135A).  In  Section  4,  the  model  is  compared  to  data 
taken  under  the  BAMM  Program.  In  Section  5,  the  results  are  summarized  and 
some  Implications  for  downward  looking  systems  are  indicated. 


2.  CLOUD  MODEL 


2.1  Viewing  Scenario 

When  viewing  the  earth  background  from  a  high  altitude  satellite  or 
balloon  platform,  the  upwelling  IR  radiation  can  be  described  in  terms  of 
two  signature  elements:  earth-atmosphere  thermal  emission  and  solar  radia¬ 
tion  scattered  into  the  instrument  line  of  sight.  The  intensity  of  the 
measured  signal  is  a  sensitive  function  of  the  sun-earth-observer  viewing 
geometry.  Clouds  complicate  the  description  of  the  radiation  field  by 
attenuating  radiation  from  the  earth  surface  and  strongly  scattering  solar 
radiation.  The  schematic  shown  in  Fig.  1  illustrates  the  viewing  scenario 
and  Indicates  the  important  IR  phenomena  which  must  be  incorporated  when 
modeling  scenes  that  contain  clouds.  The  cloud  model  described  in  these 
sections  is  designed  to  calculate  the  upwelling  IR  radiation  for  a  wide 
range  of  viewing  scenarios  which  Include  both  optically  thin  and  thick 
clouds.  When  validated  against  the  BAMM  data,  the  model  can  be  used  as  the 
basis  for  systems  design  calculations.  The  purpose  of  the  model  is  to  give 
reasonable  results  with  sufficient  flexibility  for  incorporating  different 
parameters  such  as  cloud  type  and  altitude,  atmospheric  conditions,  solar 
position,  observer  altitude,  and  background  conditions.  Since  the  amount 
of  solar  radiation  scattered  by  clouds  varies  considerably,  the  model  must 
include  the  effects  of  important  microphysical  cloud  properties  such  as 
particle  size,  distribution,  and  scattering  properties.  Additionally,  cloud 
can  be  composed  of  either  water  droplets  or  ice  particles,  or  a  mixture  of 
the  two.  The  present  model  is  designed  to  incorporate  these  features  and 
allow  the  user  to  select  from  a  range  of  cloud  properties. 
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Figure  1.  Schematic  Showing  Field  of  View  Elements  Which  Contribute 
to  the  Upwelling  Infrared  Radiation  Observed  from  the 
BAMM  Platform. 


2.2  Cloud  Model 

Considerable  theoretical  and  empirical  effort  has  been  expended  in  develop¬ 
ing  calculational  models  which  are  capable  of  predicting  the  upwelling  radia¬ 
tion  field  in  the  presence  of  clouds.  The  major  thrust  of  these  efforts  has 
been  to  develop  global  and  local  climatological  models  for  predicting  how  clouds 
affect  the  planetary  albedo,  the  surface  cooling  rate, and  the  upwelling  radiance 
data  taken  by  meteorological  satellites.  The  emphasis  of  these  studies  has  been 
on  predicting  the  net  upward  flux  (the  radiance  integrated  over  wide  bandpasses 
and  the  upper  hemisphere) .  The  Air  Force  interest  in  surveillance  problems 
dictates  that  new  effort  be  directed  toward  prediction  of  directional  radiances 
(FOV  observed  by  a  satellite  system)  within  narrow  bandpasses  (target  identifica 
tion  by  spectral  signature). 

The  calculational  models  can  be  divided  into  two  categories;  exact 

but  computationally  slow  models  and  approximate,  computationally  fast 

models.  Because  a  satellite  platform  will  view  a  large  complicated  earth 

scene  that  includes,  in  addition  to  clouds,  surface  features,  surface 

vehicles,  varying  atmospheric  conditions,  changing  sun-FOV  background- 

observer  geometry,  and  moving  targets,  the  emphasis  of  the  present  modeling 

effort  was  placed  on  the  faster  calculational  techniques.  Comparisons 

to  field  data  and  to  the  detailed  models  will  aid  in  validating  the  model 

and  quantifying  its  limitations.  Multiple  scattering  effects  within  a 

cloud  are  the  most  difficult  to  model;  several  of  the  techniques  include 

(1234) 

the  doubling  method  first  suggested  by  van  de  Hulst,  *  *  *  ’  Monte 

"'‘van  de  Hulst,  H.C.,  "A  New  Look  at  Multiple  Scattering,"  Institute  for 
Space  Studies  Report,  NASA,  New  York  (1963). 

2 

Hansen,  James  E.  and  Pollack,  James  B.,  "Near  Infrared  Scattering  by 
Terrestrial  Clouds,"  J.  Atmo.  Sci. ,  27 ,  266  (1970). 

3 

Hansen,  James  E.,  "Exact  and  Approximate  Solutions  for  Multiple  Scattering 
by  Hazy  Atmospheres,"  J.  Atmo.  Sci.,  26,  478  (1969). 

4 

Welch,  R.  M.  and  Cox,  S.  K.,  "Absorption  of  Solar  Radiation  in  Clcuds;  The 
Effect  of  Monomodal  and  Bimodal  Size  Distributions,"  Third  Conference  on 
Atmospheric  Radiation,  American  Met.  Soc.  (1978). 


Carlo  calculations,  and  multiple  stream  approaches.  The  model 

outlined  below  Is  based  on  a  two-stream  (upward  and  downward  radiance)  approach 
that  Is  a  simpler  form  of  the  multiple  stream  models.  This  approach  draws  on 
results  from  previous  efforts  in  forming  a  calculational  model  which  applies 
with  reasonable  accuracy  to  different  cloud  altitudes,  types  and  thicknesses. 

Clouds  strongly  scatter  radiation  in  the  infrared  because  of  their  high 

density  of  aerosols  (either  liquid  ^0  or  ice  crystals).  Low  altitude  cumulus 

clouds  are  usually  composed  of  liquid  1^0  or  a  liquid-ice  mixture.  (^*11) 

At  higher  altitudes  with  their  colder  temperatures,  cloud  droplets  are 

(8) 

almost  exclusively  ice.  '  Thin,  high  altitude,  cirrus  clouds  transmit 
much  of  the  upwelling  radiation  from  below,  while  scattering  significant 
amounts  of  the  impinging  solar  radiation.  This  effect  can  lead  to  enhanced  1R 
signals  in  the  shorter  wavelength  regions,  while  maintaining  information 
about  signals  arising  in  the  lower  altitudes. 


Davies,  Roger,  "The  Three-Dimensional  Transfer  of  Solar  Radiation  in  Clouds, " 
Ph.D.  Thesis,  University  of  Wisconsin,  Madison,  Wisconsin  (1976)  and 
J.  Atmo.  Sci. 

SlcKee,  T.  B.  and  Cox,  S.  K.,  "Scattering  of  Visible  Radiation  by  Finite 
Clouds,"  J,  Atmo.  Sci.,  31,  1886  (1974). 

7Liou,  Kuo-Nan  and  Wittman,  Gerald  D.,  "Parameterization  of  the  Radiative 
Properties  of  Clouds,"  J .  Atmo ,  Sci ■ ,  36,  1261  (1979). 

8 

Roewe,  Douglas  and  Liou,  Kuo-Nan,  "Influence  of  Cirrus  Clouds  on  the 
Infrared  Cooling  Rate  in  the  Troposphere  and  Lower  Stratosphere,"  J.  Appl. 
Meteor. ,  17,  92  (1978). 

9 

Harshvardhan ,  "Perturbation  of  the  Zonal  Radiation  Balance  by  a  StratospLu  r  ic 
Aerosol  Layer,"  J.  Atmo.  Sci.,  36,  1274  (1979). 

*^Pruppacher ,  Hans  R.  and  Klett,  James  D.,  "Microphysics  of  Clouds  and 
Precipitation,"  Reidel  Publishing  Co.,  Boston  (1978). 

*^Kuhn,  P.M.  and  Weickmann,  H,  K.,  "High  Altitude  Radiometric  Measurements 
of  Cirrus,"  J .  Appl .  Meteor . ,  6,  147. 
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The  upwelling  radiation  described  by  the  model  is  obtained  by  solving  the 


generalized  monochromatic  equation  of  radiation  transport 


(12,13) 


which  is: 


< 

where  the  rate  of  change  of  spectral  radiance  I(t,p)  at  an  optical  depth  t  ana 
zenith  direction  cosine  p  results  from  the  contributions  of  the  radiation  elements 
on  the  right  hand  side.  In  order,  these  terms  describe: 

•  Attenuation  of  the  radiation  due  to  absorption  and  scattering  out 
of  the  observer's  line  of  sight  (LOS); 

•  Ambient  radiation  scattered  into  the  LOS  from  all  directions; 

•  Thermal  emission  from  the  cloud  and/or  atmosphere,  and 

•  Radiation  scattered  into  the  observer’s  LOS  from  external  sources. 


The  external  source  in  the  present  problem  is  the  solar  irradiance.  The  up- 
welling  radiation  from  the  earth  and  atmosphere  is  included  as  a  boundary 
condition  for  the  lower  edge  of  a  cloud  because  it  is  a  diffuse  source  term. 
Since  the  integro-dif f erential  equation  cannot  be  solved  in  closed  form, 
approximate  solutions  that  describe  the  salient  features  must  be  developed. 
Before  describing  these  approximate  solutions  to  the  radiative  transfer 
equation,  model  approximations  to  the  process  of  scattering  by  single  particles 
are  described  and  developed. 


^Chandrasekhar ,  Radiative  Transfer,  Dover,  New  York  (1950), 


JMc Elroy ,  M.B.,  "The  Composition  of  Planetary  Atmospheres," 
Spectrosc.  Radiat.  Transfer.  11,  813  (1971). 
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2.2.1  Single  Scattering 


The  single  scattering  event  is  described  by  the  differential  cross' 
section  for  an  Incident  photon  to  be  scattered  into  an  arbitrary  direction. 
Thi6  cross-section  is  written  as  a  phase  function,  P(yg),  which  is  defined  as 


*<*>.> 


A  d° 

a  dy 

s 


(2) 


where  is  the  differential  cross  section  for  a  given  particle  size  dis¬ 
tribution,  and  yg  is  the  cosine  of  the  scattering  angle.  Scattering  calcula¬ 
tions  based  on  Mie  theory  provide  a  good  description  of  the  scattering  phase 
function  for  spherical  particles.  ^ The  first  step  in  these  calcu¬ 
lations  is  to  determine  the  distribution  of  droplet  sizes.  Figures  2  and  3 

give  examples  of  the  scattering  phase  function  (calculated  from  Mie  theory) 

(2) 

for  different  particle  size  distributions.  Figure  2  shows  the  single 

scatter  phase  function  for  moderate  and  large  mean  particle  sizes.  The 

curves  have  been  displaced  vertically  by  factors  of  ten.  The  horizontal 

lines  on  the  curves  correspond  to  the  point  P(0)  »  1.  Figure  3  shows  the 

single  scatter  albedo  for  spherical  water  and  ice  particles  as  reported  by 

(2) 

Hansen  and  Pollack.  Two  different  mean  particle  radii  are  shown,  2  ym  and 
32  ym.  The  particle  size  distributions  are  given  by^^ 


n(r)  «  r 


6 


exp  (-6r/r  ) 
m 


(3) 


14 

Friedlander,  S.K.,  Smoke,  Dust  and  Haze,  John  Wiley  and  Sons,  New  York 

(1977). 

^Zachor,  A.S.,  Holzer,  J.H.,  Smith,  F.G.,  "IR  Signatures  Study,"  Honeywell 
Electro-Optics  Center  Report  No.  7812-8,  December  1978. 

^Shettle,  E.  (Calculations  reported  in  Zachor,  et  al.,  Ref.  15). 

^Deirmendjian,  "Scattering  and  Polarization  Properties  of  Water  Clouds  and 
Hazes  in  the  Visible  and  Infrared,"  Appl.  Opt.,  187  (1964). 
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e  2.  Single-Scattering  Phase  Function  for  Two  Size  Distributions 
of  Spherical  Particles  Showing  the  Effect  of  Absorption. 

The  Horizontal  Lines  Correspond  to  the  Value  1.0.  From 
Ref.  (2). 
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Figure  3.  Comparison  of  the  Single  Scattering  Albedo  For  Water 
Droplets  and  Spherical  Ice  Crystals  at  Two  Different 
Mean  Particle  Radii.  From  Ref.  (2). 


where  r  is  the  mean  particle  radius.  Since  conventional  Mie  calculations 

ID 

strictly  apply  only  for  spherical  particles,  a  basic  assumption  of  virtually 
all  cloud  models  is  that  the  particles  are  spherical.  Hie  type  calculations 

have  been  carried  out  for  cylindrical  particles,  which  are  more  represents- 

(18)  (18) 
tive  of  ice  crystals.  '  Figure  4,  which  is  taken  from  Liou's  work, 

indicates  that  the  cylinders  scatter  a  little  more  than  the  spheres  and  have 

less  structure  at  larger  scattering  angles.  The  assumption  of  spherical 

particles  is  very  good  for  water  droplets  but  is  less  valid  for  ice  crystals 

which  are  large  (10-100  ym)  and  often  hexagonal  in  shape.  From  these  curves 

one  sees  that  the  increase  in  the  back  scatter  becomes  weaker  as  the  particle 

shapes  become  more  asymmetric.  Therefore,  we  have  not  made  any  attempt  to 

include  a  back-scattered  enhancement  in  the  present  model  because  higher 

altitude  clouds  contain  ice  crystals  with  their  nonspherical  shapes. 


Figure  4.  Scattering  Phase  Function  for  Uniformly  Oriented  Ice  Spheres  at 
Three  Different  Wavelengths.  The  Curves  are  Displaced  for 
Clarity  and  the  Horizontal  Lines  Through  the  Curves  Correspond 
_ _ to  the  Value  1.0.  From  Ref.  (18). 

18 

Liou,  Kuo-Nan,  "Light  Scattering  by  Ice  Clouds  in  the  Visible  and  Infrared 
A  Theoretical  Study,"  J.  Atmo.  Sci..  29,  524  (1972). 


Detailed  Mie  calculations  for  every  particle  else,  wavelength,  and 
scattering  angle  are  time  consuming  and  beyond  the  scope  of  an  engineering 
type  model.  Therefore,  a  function  which  approximates  the  phase  functions 
derived  from  Mie  calculations  is  used  in  the  model  development.  The  most 
widely  used  and  studied  phase  function  is  the  Henyey-Greenstein  (H-G) 
function,  which  was  first  Introduced  in  1941 


PH-G(C0S  V 


(1  +  g2  -  2g  cos  e)3/2 

s 


where  cos  0g  is  the  cosine  of  the  scattering  angle.  The  advantage  of 

this  scattering  function  over  detailed  results  from  Mie  theory  is  that 

it  provides  a  reasonable  one-parameter  fit  to  the  detailed  calculations  for 

the  angular  distribution  of  light  scattering  by  spherical  droplets.  The 

parameter  is  the  asymmetry  factor  g,  which  is  defined  as  the  mean  value  of 

Vi  ,  the  cosine  of  the  scattering  angle.  Figure  5  shows  a  comparison  of  the 
8 

Henyey-Greenstein  representation  with  Mie  theory.  The  H-G  curve  gives  a 

reasonable  overall  description,  although  it  does  not  have  the  small  back 

scattering  increase  of  the  Mie  calculation.  The  back  scattering  peak  could 

be  easily  parameterized  with  a  second  H-G  phase  function.  Figure  6,  which 

(31 

is  also  taken  from  the  work  of  Hansen,  shows  a  similar  comparison  for  exact 
calculations  of  cloud  directional  radiance  for  two  different  optical  depths 
and  solar  angles.  In  the  present  model,  spectral  variation  in  the  phase 
function  is  Introduced  by  varying  the  asymmetry  parameter  g,  based  on  published 
Mie  calculations 

Molecular  absorption  by  the  atmospheric  gases  within  clouds  moderates 
the  Intensity  of  multiple  scattering  within  the  clouds  for  those  spectral 
bands  where  the  gases  absorb  strongly.  In  these  spectral  regions  a  photon 

T9  ~ 

Henyey,  L.G.  and  Greenstein,  J.L.,  "Diffuse  Radiation  in  the  Galaxy," 

Aatrophys.  J..  93,  70  (1941). 
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Relative  Intensity 


Figure  5.  Comparison  of  the  Angular  Distribution  of  the 
Henyey-Greenstein  Phase  Function  to  Calcula¬ 
tions  from  Mie  Theory.  From  Ref.  (3). 


Figure  6.  Comparison  of  the  Angular  Distribution  at  Two  Optical 
Depths  from  A  Cloud  For  Two  Solar  Azimuth  Angles.  An 
Exact  Calculation  with  Mie  Theory  is  Compared  to  One 
With  the  Henyey-Greenstein  Phase  Function.  The  Wave¬ 
length  is  0.7  um.  From  Ref.  (3). 


will  have  little  opportunity  for  more  than  one  or  two  scattering  events  before 
being  absorbed.  The  effect  of  molecular  absorption  is  handled  by  modifying  the 
single  particle  scattering  albedo,  which  is  defined  as  the  ratio  of  the  scatter 
ing  to  extinction  cross  sections.  Under  this  assumption  the  albedo  u  is  given 
by 


o 


where  Os  is  the  single  particle  volume  scattering  coefficient,  is  the 
particle  extinction,  and  is  the  molecular  volume  absorption  coefficient. 
Thus  the  effective  single-particle  scattering  is  reduced  as  the  molecular 
absorption  coefficient  becomes  larger. 

2.2.2  Multiple  Scattering 

The  multiple  scattering  term,  which  is  the  second  one  in  the  radiation 
transport  equation,  Eq.  (1),  describes  the  contribution  of  ambient  radiation 
scattered  into  the  observer's  LOS.  It  is  this  term  which  makes  Eq.  (1) 
unsolvable  in  closed  form.  Highly  accurate  numerical  solutions  can  be 
obtained  using  Monte  Carlo  calculational  techniques.  Such  calculations, 
which  follow  the  path  of  every  photon,  are  not  suitable  for  a  quick, 
flexible  model  with  its  emphasis  on  describing  the  overall  FOV  scene  with 
allowances  for  variations  in  the  scene  and  viewing  parameters.  However,  Monte 
Carlo  calculations  are  extremely  useful  as  a  check  on  more  approximate 
models. 

The  model  which  we  selected  approximates  the  multiple  scattering  integral 
(13  18) 

by  two  terms;  ’  '  These  terms  can  be  considered  as  the  upwelling 

and  downwelling  radiance.  The  Integral  o’er  y'  or  cos  6'  in  Eq.  (1)  is 
approximated  by  integrals  over  the  upper  and  lover  he^isnberes.  with  the 
assumption  of  isotropic  intensity  in  each  hemisphere.  Once  this  assumption 
is  made,  Eq.  (1)  becomes  a  set  of  two  coupled  differential  equations  for  the 
upward  and  downward  radiance. 


The  two-stream  modeling  approach  has  been  used  extensively  in  the  cloud 

literature,  especially  in  Refs.  (7)  to  (10),  and  variations  in  detailed 

approaches  are  reviewed  in  Ref.  (20).  These  studies  are  best  summar iatJ 

by  saying  that  the  two-stream  model  gives  good  results  for  the  total 

upwelling  radiance  and  adequate  results  for  directional  scattering.  Although 

detailed  comparisons  with  exact  scattering  models  have  not  been  made  for 

arbitrary  scattering  geometries,  we  have  found  that  the  two-stream  approach 

predicts  directional  results  which  are  in  good  agreement  with  an  exact  Monte 
(21) 

Carlo  model,  for  the  case  of  zenith  sun. 

2.3  Two-Stream  Model 

The  formulation  of  the  present  cloud  model  is  based  on  several  assump¬ 
tions.  A  cloud  of  arbitrary  and  varying  optical  thickness  is  modeled  by  a 
plane  parallel  atmosphere  in  which  the  single  scattering  albedo  (u>v)  is 
Independent  of  depth.  Local  thermodynamic  equilibrium  between  particles 
and  atmospheric  gases  is  assumed,  with  a  uniform  temperature  throughout  the 
cloud.  With  these  assumptions,  the  equation  of  radiative  transfer  (with 
total  extinction  optical  depth  T  as  the  independent  variable)  is 

1  2TT 

XV  I  (T.y’.r)  dd>’  dll’ 

-l  o 

-  (l  -  u>v)  bv  (T)  -  exp  (-t/yo)  P(y,4>;  -  ,4>o)  0* 


Meador,  W.E.,  and  Weaver,  W.R.,  "Two-Stream  Approximations  to  Radiative 
Transfer  in  Planetary  Atmospheres:  A  Unified  Description  of  Existing 
Methods  and  a  New  Improvement,"  J.  Atmo.  Sci.  37,  630  (1980). 
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Bernstein,  L.S.,  Duff,  J.S.,  Stanton,  A.C.,  and  Robertson,  D.C.,  in 

preparation. 


this  equation  may  be  rewritten  in  terms  of  a  source  frunction  J  (t)  asv  ’  ' 


I  -  J 
v  v 


Equation  (5)  may  be  solved  formally  for  the  emergent  radiance  at  the  top  of 
the  cloud.  The  formal  solution  is 


Iv  (t  -  0,y) 


-T./u 

Wy)e  + 


f1  e~l/v 

J  y  dT, 


where  Tj  is  the  optical  thickness  of  the  cloud.  Because  the  source  function 
includes  the  multiple  scattering  term  which  involves  the  unknown  radiance, 
an  approximate  method  is  adopted  for  solving  the  radiative  transfer  equation. 
This  method  consists  of  two  steps.  In  the  first  step,  the  two-stream  approx¬ 
imation  is  made  and  the  multiple  scattering  integral  is  represented  by  integrals 
over  the  upper  and  lower  hemispheres.  In  the  second  step,  the  approximate 
multiple  scattering  solution  obtained  in  the  first  step  is  combined  with  the 
thermal  emission  term  and  an  exact  representation  of  the  single  scattering 
term.  These  terms  form  the  source  function,  which  is  used  in  the  formal 
solution,  Eq.  (6). 

The  azimuthally  averaged  form  of  the  radiative  transfer  equation  is  used 
in  developing  the  approximation  for  multiple  scattering.  This  form  of  the 
equation  is 


dl  _  u> 

U  T-  *  1  - 


(1  -  w)B  -  —  exp(-i/u  )  P(p,-  u  ) 


P(y,  p')  I(vi')  dv.  ’ 


vhere 


P(cos  0)  d(cos  0) 


1. 


-(8) 


The  two-stream  approximation  is  made  by  assuming  that  the  radiance  is 
Independent  of  zenith  angle  in  each  of  the  upper  and  lower  hemispheres,  i.e 

I(u)  =  I+,  0  <  p  <  1 

I(p)  =  I“,  -1  <  u  <  0  .  (9) 

Equation  (7)  may  then  be  integrated  over  the  upper  and  lower  hemispheres  to 
form  two  coupled  equations  in  1+  and  1“ : 


1 

2  dT 


I+  -  |  <1+  +  I~)  -  ~  (1+  -  I")  -  (1  -  io)B 


^  exp(-t/pQ)  u[l  -  G(uo)  J  , 


1  d7”  *  I'  -  f  (I+  +  I")  +  f  -  (I+  ‘  1_>  ‘  O  "  u)B 


-  j  exp(-T /uo)  u>  £  1  +  G(pq)  J  , 


where  G 


-  I  P(u.  p')  dp'  -  1  , 

»  r\ 


H-f 


G(p)  dp 


The  quantities  G(p)  and  H,  which  involve  hemispherical  integrals  of  the 
scattering  phase  function,  have  been  evaluated  for  the  Henyey-Greenstein 
function  in  Ref .  (22) . 


The  boundary  conditions  for  the  problem  are  given  by  continuity  of  the 
radiance  at  the  cloud  top  (r  15  0)  and  cloud  bottom  (t  =  tj): 

1~(t  *  0)  =  0  ,  (14a) 

l+(1  -  T:)  =  1B(W)  .  (14b) 

Equation  (14a)  states  that  the  diffuse  downward-directed  radiance  at  the 

cloud  top  is  zero.  Equation  (14b)  equates  the  upward-directed  diffuse 

radiance  at  the  cloud  bottom  to  I  (p),  which  represents  the  upwelling  radiance 

B 

from  the  ground  (thermal  emission  and  solar  scatter)  plus  atmospheric  emission 
along  the  path  from  the  ground  to  the  cloud. 

Equations  (10)  and  (11)  with  boundary  conditions  (14a)  and  (14b)  are 
sufficient  to  solve  for  1+  and  1  ,  the  initial  approximations  to  the  upward 
and  downward  diffuse  radiance.  The  details  of  this  solution  are  algebraically 
tedious  and  will  not  be  3iven  here.  After  this  solution  is  obtained,  the 
approximate  source  function  is  formed  as 

JV(T)  =  |  (I+  +  1")  +  |  G(u)(I+  -  I")  +  (1  -CJ)BV(T) 

+  ^  w  exp(--r/uo)  PH_G<U.  <fQ)  •  (15) 

With  an  approximation  to  the  source  function  now  defined,  the  solution 
for  the  emergent  radiance  at  the  top  of  the  cloud  mav  be  found  from  Eq .  (6). 
This  solution  is  given  here  for  completeness: 

_ 

Wiscombe,  W.J.,  and  Grams,  G.W.,  "The  Backscattered  Fraction  in  Two-Strear 
Approximations,"  J.  Atmo,  Sci.  33,  2440  (1976). 
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In  this  equation,  u  is  the  cosine  of  the  observer  zenith  angle,  u  is  the 
cosine  of  the  solar  zenith  angle,  Tj  is  the  total  extinction  optical  thickness 
of  the  cloud  in  the  vertical  direction,  ttF  is  the  solar  irradiance  at  the  tcp 
of  the  cloud,  B  (T)  is  the  Planck  function,  and  I0  is  the  upwelling  radiance 
at  the  lower  boundary  of  the  cloud.  Other  terms  in  this  equation  are: 
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2-4  Implementation  of  the  Cloud  Model 

The  cloud  model  is  implemented  as  a  modification  to  the  AFGL  atmospheric 
transmittance  code  L0WTRAN3B.'^^  A  cloud  layer  of  arbitrary  altitude  and  spatial 
thickness  is  defined  in  the  LOWTRAN  input  atmospheric  profile  by  specifying 


Selby,  J.E.A.,  Shettle,  E.P.,  and  McClatchey,  R.A.,  "Atmospheric  Transmit¬ 
tance  from  0.25  to  28.5  ym:  Supplement  LOWTRAN 3B,"  AFGL-TR-76-0258,  Air 
Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731  (November  1976). 
ADA040701, 
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the  aerosol  number  density,  temperature,  and  water  vapor  concentration  In  the 
cloud.  Standard  model  atmospheres  or  radiosonde  data  are  used  to  specify- 
the  rest  of  the  atmospheric  profile.  Aerosol  extinction  and  absorption 
coefficients  which  pertain  to  liquid  water  droplets  or  ice  particles  are 
used  instead  of  the  properties  corresponding  to  the  standard  LOWTRAN  atmos¬ 
pheric  aerosol  models.  These  wavelength-dependent  coefficients  are  taken 
from  Mi e  theory  calculations  by  Shettle,^*^  as  reported  by  Zachor,  et  al.^1^ 
A  log-normal  particle  size  distribution, 


dN(r) 

dr 


(log  r  -  log 
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rm)2 
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was  assumed  in  Shettle's  calculations,  and  results  are  available  for  both  water 
and  ice  for  a  range  of  mean  particle  si2es.  The  modified  LOWTRAN  code  combines 
these  aerosol  parameters  with  its  molecular  absorption  band  model  to  calculate 
the  transmittance  through  the  cloud  layer. 

A  special  subroutine  is  added  to  the  code  to  calculate  the  cloud  radiance, 
including  cloud  thermal  emission,  scattered  solar  radiation,  and  scattered 
or  directly  transmitted  upwelling  radiance  from  below  the  cloud.  This  cal¬ 
culation  is  an  evaluation  of  Eq.  (16)  above.  The  calculation  of  cloud 
radiance  from  Eq.  (16)  requires  values  for  the  wavelength  dependent  quantities 
wv»  8.  F.  T1»  and  1^.  The  single  scattering  albedo  is  calculated  as 
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where  np°gcat  and  np°ext  are  the  volume  particle  scattering  and  extinction 
coefficients  and  n  0  is  the  volume  gas  absorption  coefficient.  Molecular 

O  © 

absorption  effects  included  in  the  cloud  model  are  therefore  characterized 
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by  the  LOWTRAN  transmittance  model,  i.e*t  band  and  continuum  absorption  by 
Important  atmospheric  gases  are  considered,  with  a  spectral  resolution  of 
20  cm  The  anisotropy  parameter,  g,  which  characterizes  the  single  scat¬ 
tering  phase  function,  is  specified  as  input  to  the  model  using  the  results 
of  the  Mie  theory  calculations  by  Shettle/16^  The  solar  flux  ttF  incident 
at  the  cloud  top  is  given  by 
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where  Hq(v)  is  the  solar  lrradiance  at  the  top  of  the  atmosphere  and 
Tv (°°  ■*  cloud)  is  the  atmospheric  transmittance  from  the  top  of  the  atmosphere 
to  the  cloud  top,  calculated  by  LOWTRAN.  The  cloud  optical  thickness,  T^(v), 
may  be  obtained  from  the  LOWTRAN  transmittance  calculations  and  includes 
the  effects  of  gaseous  absorption,  aerosol  absorption,  and  aerosol 
scattering.  The  diffuse  spectral  radiance  at  the  bottom  of  the  cloud, 

IQ,  is  calculated  by  LOWTRAN  with  the  addition  of  solar  scatter,  and  is 
given  by 


I0(v) 


•  [■ 

*  I 


cos0. 


e  B  (T  ) 
g  v  g 


r  H_ 
g  0 


>]• 


T  (“  -*•  ground)!  •  T(ground  -*•  cloud) 


cloud 


V”  fe 


(26) 


ground 


Equation  (26)  includes  the  effects  of  ground  thermal  emission,  sunlight 
scattered  from  the  ground,  and  emission  from  atmospheric  gases  beneath  the 


cloud.  The  ground  emissivity  e  (v)  and  reflectivity  r  (v)  are  input 

g  g 


parameters  in  the  model. 


The  apparent  cloud  signal  at  observer  altitude  is  calculated  as 

rbserver 

BV(T)  dz  .  (2 

oud 

In  this  expression,  Iclouj(v)  is  the  cloud  source  radiance  calculated  accord¬ 
ing  to  Eq.  (16)  as  described  above,  t(cloud  -►  observer)  is  the  atmospheric 
transmittance  between  the  cloud  top  and  the  observer,  and  the  integral  term 
is  the  emission  from  atmospheric  gases  between  the  cloud  and  the  observer. 

To  summarize  this  section,  the  cloud  model  is  combined  with  the  LOWTRAN 
atmospheric  transmission  model  to  form  a  general  atmospheric  radiative 
transfer  code  which  is  applicable  to  the  calculation  of  infrared  background 
signals  (Including  clouds)  which  could  be  detected  by  high  altitude  downward¬ 
looking  sensors.  The  inputs  to  the  model  are  background  characteristics 
(surface  altitude,  temperature,  emisslvity,  and  reflectivity),  cloud  charac¬ 
teristics  (altitude,  thickness,  temperature,  water  vapor  concentration, 
particle  type,  particle  number  density,  and  mean  particle  size),  observer 
altitude,  and  observer  and  solar  angles  (zenith  and  azimuth).  Standard 
atmospheric  models  or  radiosonde  data  may  be  used  to  define  the  atmospheric 
profiles  of  temperature,  water  vapor  concentration,  etc.  The  output  from 
the  model  is  the  calculated  spectral  radiance  (20  cm  *  resolution)  at  the 
observer's  position  due  to  thermal  emission  from  clouds  and  the  earth  surface, 
sunlight  scattered  by  clouds  and  earth  Burface,  and  thermal  emission  from 
the  atmosphere.  Integrated  signals  may  be  calculated  for  specified  bandpasses 


1: 5  Other  Cloud  Modeling  Work 


Two  previous  cloud  modeling  efforts,  which  have  directly  addressed  the 

prediction  of  directional  radiances  in  spectral  bandpasses  of  possible  interest 

(24) 

for  surveillance  applications,  are  a  study  by  S.J.  Young  at  the  Aerospace 
Corporation,  and  cloud  modeling  work  at  Honeywell  by  Zachor,  Holzer,  and 
Smithf*^  Young's  work  ^  reviews  several  approximate  schemes  for  calculating 
scattered  radiation  from  clouds  and  adopts  a  method  which  is  an  empirical 
modification  of  a  single-scattering  model.  Comparisons  with  exact  solutions 
are  given.  In  the  work  of  Zachor,  et  al.f*^  an  approximate  two-stream 
solution  to  the  radiative  transfer  equation  is  developed,  and  the  solution 
is  compared  with  the  cloud  data  obtained  in  the  2.5  -  6  vm  region  by 
the  AFGL  KC135A  Flying  Laboratory.  The  model  discussed  in  this  report  is 
similar  to  the  model  of  Zachor,  et  al.,  with  the  important  difference  that 
molecular  absorption  by  atmospheric  gases  is  included  in  the  model.  In  addition, 
the  present  model  ensures  that  the  proper  amount  of  energy  is  scattered  into  the 
upper  and  lower  hemispheres.  Incorporating  the  present  model  in  LOWTRAN  permits 
an  evaluation  of  atmospheric  attenuation  and  background  source  effects  on  the 
apparent  cloud  signature.  These  effects  are  very  important  in  calculating 
cloud  radiance  within  molecular  absorption  bands  and  for  predicting  the  effects 
of  cloud  altitude  on  apparent  cloud  signal. 
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3.  COMPARISONS  WITH  CLOUD  SPECTRAL  DATA 


In  this  section  the  cloud  radiance  model  is  compared  to  cloud  IR  sig¬ 
nature  data  taken  from  two  different  aircraft  viewing  platforms.  The 
observations  were  made  in  the  downward  looking  mode  from  above  a  solid 
cirrus  cloud  deck  in  the  2. 5-6.0  ym  spectral  region.  The  measurement 
platforms  were  a  U2  aircraft  at  15  km  altitude  and  a  KC135A  aircraft  at 
5-10  km.  The  IR  signature  in  the  shorter  wavelength  region  is  dominated  by 
solar  radiation  singly  and  multiply  scattered  by  the  cloud  deck.  At  longer 
wavelengths  thermal  emission  from  the  cloud  Itself  becomes  increasingly 

Important.  The  next  subsection  discusses  comparisons  of  model  calculations 

(25) 

for  data  taken  by  Arthur  D.  Little  Co.  from  a  U2  aircraft,  and  the  second 
subsection  compares  model  calculations  to  AFGL  data  taken  from  a  KC135A 
aircraft  Both  of  these  measurements  have  the  same  general  viewing 

geometry,  a  schematic  of  which  is  shown  in  Fig.  7.  The  sun  was  near  the 
horizon  so  these  data  are  for  small  total  scattering  angles  (6g  -  10-15°, 
as  defined  in  Eq.  (A)).  One  data  comparison  is  made  with  the  incident 
solar  direction  and  measurement  direction  at  right  angles  so  that  0g  =  89°. 

3.1  2.7  ym  H20/C02  Band 

Measurements  were  taken  over  a  cumulus /cirrus  cloud  deck  with  its  top 

(25) 

at  an  altitude  of  15.25  km.  The  measurements  were  made  from  a  U2  flying 
300  m  above  the  clouds.  A  spectroradiometer  was  used  to  take  data  in  the 
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Espinola,  R.P.,  "Spectral  and  Spatial  Properties  of  Cloud  Backgrounds 
From  2.65  to  2.95  Microns,"  Report  No.  70505F,  Arthur  D.  Little,  Inc., 
Cambridge,  MA  02140  (1970). 
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Figure  7.  Schematic  Showing  the  Measurement  Geometry  for  Radiation 
From  the  Top  of  A  Uniform  Horizontal  Cloud  Deck.  In 
This  Example  the  Observing  Aircraft  is  Looking  Towards 
the  Sun,  so  the  Aspect  Angle  <|>  Between  the  Incident 
Solar  Direction  and  Exiting  Radiation  is  Zero. 

2.65-2.95  um  region,  with  a  spectral  bandpass  of  0.02  um.  The  uncertainty 

in  this  spectral  region  is  10%  for  each  measurement  with  a  minimum  uncertainty 
—6  2 

of  approximately  10  watts/sr/cm  /pm.  At  the  time  of  measurement  the  sun 
was  near  the  horizon  with  a  zenith  angle  of  88°.  The  U2  aircraft  was  looking 
down  from  an  altitude  of  15.5  km  or  approximately  1000  ft  above  the  cloud 
deck  with  an  observer  zenith  angle  of  80°.  The  aircraft  was  looking  almost 
directly  towards  the  sun  with  a  relative  polar  angle  of  1°.  The  scatter¬ 
ing  angle  between  the  incident  and  outgoing  radiation  is  12°. 

Figure  8  shows  a  comparison  of  data  taken  by  ADL  with  calculations 
using  the  Aerodyne  cloud  scattering  model.  The  key  features  in  the  data  are  the 
strong  COj  absorption  bands  at  2.69  and  2.77  um,  which  lead  to  high  altitude 
absorption  of  the  incident  solar  energy.  Although  there  is  no  information 
on  the  particle  type  and  size  distribution,  ice  crystals  are  certainly 
dominant  because  of  the  high  altitude.  The  calculations  assumed  a  mean 


Wavelength, 


Figure  8.  Comparison  of  Model  Calculations  to  Arthur  D.  Little  Data 
For  Top  Viewing  of  A  Cloud  in  the  2.7  um  Band.  Two  Cal¬ 
culations  with  Mean  Particle  Radii  of  4  and  32  um  are 
shown . 

(18 

size  of  32  pm  which  is  consistent  with  Liou's  cylindrical  crystal  model 

(2) 

and  with  Hansen  and  Pollack.  The  calculations  capture  the  important 
spectral  features  of  the  data,  i.e.,  CO^  absorption  and  light  scattering 
by  the  clouds,  although  they  are  consistently  a  factor  of  3-10  below  the 
data.  The  data  vary  by  a  factor  of  two,  with  Fig.  8  being  the  highest 

/  A  P  \ 

measured  values.  The  scattering  angle  of  12  is  within  the  scattering 

peak,  and  the  results  are  extremely  sensitive  to  assumptions  about  the 
single-scatter  phase  function  and  details  of  the  multiple  scattering  model 


As  discussed  In  Section  2,  the  Henyey-Greensteln  function  gives  good  overall 
predictions  for  the  intensities;  however,  it  can  be  off  by  as  much  as  a  factor 
of  2  or  more  in  its  detailed  angular  dependence  in  the  forward  scattering  peak. 
Uncertainties  in  the  phase  function  for  near-forward  scattering  are  discussed 
below  in  subsection  3.2.  Direct  multiple  scattering  contributions  (second  and 
third  scatterings)  in  the  forward  peak  region  will  lead  to  increased  radiation. 
The  model  in  its  present  formulation  includes  multiple  scattering  in  a  diffu¬ 
sion  approximation  and  direct  single  scattering.  The  calculations  presented 
here  are  based  on  the  log  normal  distributions  and  Mie  theory  calculations  for 
the  total  scattering  cross  section.  Other  calculations  with  different  par¬ 

ticle  size  distributions  (Dermendjian  Model,  with  a  mean  radius  of  4  pm) 
predict  that  the  shorter  wavelength  contributions  are  enhanced  by  10-30%.  The 
comparisons  of  phase  functions  discussed  below  show  that  a  factor  of  four  varia¬ 
tion  in  the  detailed  angular  dependence  can  occur  in  the  forward  scattering  peak. 

3.2  Comparisons  with  AFGL  Data 

Interferometric  measurements  (1  cm  *  resolution)  of  cloud  spectral  signa- 

(Ofs) 

tures  were  taken  as  part  of  the  AFGL  Aircraft  Measurement  Program,  The 

data  discussed  here  pertain  to  a  thick  uniform  cloud  deck  at  an  altitude  of 

9.3  km.  The  solar  zenith  angle  was  82°.  The  measurements  were  taken  from  an 
altitude  0.7  km  above  the  cloud,  with  an  observer  zenith  angle  of  87.5°.  The 
relative  azimuth  angles  (measured  from  the  direction  of  forward  scattering) 
are  0°  (scattering  angle  =  10°)  and  90°  (scattering  angle  =  90°).  Local 
radiosonde  data  were  used  to  characterize  the  atmospheric  conditions. 

Comparisons  of  cloud  model  calculations  with  the  AFGL  spectral  data  are 
shown  in  Figs.  9  and  10.  The  cloud  was  assumed  to  consist  of  ice  particles 
with  a  mean  diameter  of  32  pm.  The  agreement  between  the  model  and  data  is 
good  at  the  longer  wavelengths,  where  the  cloud  thermal  emission  dominates 

solar  scattering  effects.  The  strong  atmospheric  absorption  by  C0o  in  the  2200  - 

-1  1 
2400  cm  region  is  evident  in  both  figures.  At  shorter  wavelengths  where  the 
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scattered  solar  radiation  is  most  important,  the  cloud  model  overpredicts 
the  measured  signal  by  a  factor  of  2  -  4  for  the  case  of  near-forward 
scattering  (Fig.  9).  In  Fig.  10,  the  measured  and  predicted  cloud  signals 
at  shorter  wavelengths  are  lower  by  an  order  of  magnitude  than  the  forward 
scattering  case,  because  the  scattering  angle  is  >90°.  The  model  and  the 
data  are  in  good  agreement  for  this  non-forward  scattering  case.  Calculations 
performed  by  Zachor,  et  alf*^  are  also  shown  in  Figs.  9  and  10.  The  Zachor 
model  is  similar  to  the  present  model,  with  the  exception  that  the  present 
model  accounts  for  molecular  absorption  between  scattering  events  within  the 
cloud.  This  feature  is  evident  in  Fig.  9,  where  water  vapor  absorption  from 
3100  -  3500  cm  ^  leads  to  reduced  solar  scattering  in  the  predictions  of  the 
present  model. 

For  the  case  of  near-forward  scattering,  the  solar  scattering  contribution 
to  the  cloud  signal  is  dominated  by  single  scattering.  In  the  single  scattering 
limit  the  solution  to  the  radiative  transfer  equation  (optically  thick  cloud) 
reduces  to  the  result. 


I(u,  <*>) 


V  +  P 


P11-G(CCS  6s) 
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The  cloud  radiance  in  this  limit  is  thus  directly  proportional  to  the  single 
scattering  phase  function.  In  the  forward  scattering  peak,  the  phase  function 
Is  sensitive  to  several  detailed  assumptions  in  the  model  (for  example,  the 
mean  particle  size  or,  in  the  case  of  ice  crystals,  the  particle  shape).  The 
uncertainties  in  these  assumptions  give  rise  to  an  overall  uncertainty  in  the 
model  predictions,  which  represents  an  inherent  limit  in  the  accuracy  of 
idealized  models  in  predicting  real  cloud  properties. 

A  range  of  values  for  the  single  scattering  phase  function  P  £cos(l  P°)J 
appropriate  to  the  description  of  scattering  by  ice  particles  is  given  in 
Table  1.  The  table  shows  a  comparison  of  the  phase  function  from  Mit  theory 
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Figure  9.  Comparison  of  Aerodyne  and  Honeywell  Calculations  to  AFGL 

Data  for  Top  Viewing  of  a  Cloud  Deck  (Scattering  Angle  =*10°). 
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Figure  10.  Comparison  of  AFGL  Data  to  Honeywell  and  Aerodyne  Calculation 

For  Viewing  the  Same  Cloud  Deck  as  Fig.  9  (Scattering  Angle  5  ) 
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Table  1.  Comparison  of  Single  Scatter  Phase 
Functions  in  the  Forward  Scattering 
Peak  (Scattering  Angle  *  10°) 


X(pra) 

P(cos  0) 

Henyey-Greenstein 

Mie 

(27) 

Theory v; 

Shettled^) 

Liou<27) 

spheres 

cyl inders 

2.5 

1.09 

4.52 

- 

— 

3.0 

1.94 

2.10 

-2 

~8 

calculations  of  Liou^27^  for  ice  cylinders  and  ice  sphere  (30  urn  radius) 
with  the  Henyey-Greenstein  function  using  asymmetry  parameters  calculated 
by  Shettle^*^  (32  pm  radius  particles)  and  by  Liouf27^  The  comparison 
of  the  phase  functions  at  X  *  3  pm  indicates  a  factor  of  4  variation.  In 
the  absence  of  detailed  microphysical  data  on  particle  number  density,  size 
distribution,  and  particle  shape,  there  is  not  a  priori  reason  for  preferring 
any  of  the  phase  functions  in  Table  1.  A  much  broader  base  of  data  compari¬ 
sons  is  required  in  order  to  identify  the  appropriate  microphysical  data  for 
providing  a  consistent  description  of  scattering  from  clouds. 

The  alternative  phase  functions  listed  in  Table  1  would  each  result 
in  a  higher  predicted  cloi  ’  signal  for  the  case  presented  in  Fig.  9  than 
the  present  model  predictions,  which  are  based  on  the  Henyey-Greenstein 
function  using  the  Shettle^1^  asymmetry  parameters.  This  trend  would  result 
in  poorer  agreement  with  the  data  presented  in  Fig.  9.  On  the  other  hand,  the 
difference  in  the  ice  cylinder  and  sphere  phase  functions  at  >  =3.9  pm  accounts 
for  the  difference  between  data  and  calculations  for  the  high  altitude  clouds  of 
Fig.  8.  Other  effects  which  are  not  considered  in  the  model  may  be  important, 
however.  For  example,  surface  "roughness"  at  the  cloud  tor-  vhich  is  assumed 
smooth  in  the  idealized  model)  could  result  in  a  reduced  sin*-,.  scattering 
contribution  from  the  top  layers  of  the  cloud  for  the  grazing  geometry 

27Liou,  K.N.,  "Transfer  of  Solar  Irradiance  through  Cirrus  Cloud  Lavers,' 

J.  Geophys.  Res.,  78 ,  1A09  (1973). 
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applicable  to  Fig.  9.  As  Indicated  in  Fig.  11,  surface  roughness  would  result  in 
attenuation  of  much  of  the  forward  scattered  radiation,  resulting  in  a  signifi¬ 
cantly  smaller  cloud  signal.  A  detailed  consideration  of  these  effects  was 
outside  the  scope  of  the  present  study.  For  purposes  of  the  present  work,  the 
two-stream  cloud  model  shows  good  qualitative  agreement  with  data  in  the  spectral 
dependence  of  cloud  signals.  For  near  forward  scattering  geometries,  however, 
quantitative  agreement  is  sensitive  to  the  detailed  cloud  microphysics  and 
local  scattering  geometry. 


Aiao-izo 


Figure  11.  Attenuation  of  Singly  Scattered  Radiation  by  Cloud  Surface  Roughness. 
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4.  COMPARISON  OF  MODEL  PREDICTIONS  WITH  BAMM  DATA 


In  this  section,  the  Aerodyne  cloud  and  background  radiance  model 

discussed  in  the  previous  sections  is  applied  to  the  analysis  of  field  data 

(28) 

taken  in  the  BAMM  program.  '  These  data  consist  of  radiometer  measurements  in 
the  2. 7-3.1  pm  region  from  the  BAMM  test  flight  of  September  28,  1978  at 
Holloman  Air  Force  Base,  New  Mexico.  Model  predictions  are  compared  with 
BAMM  measurements  for  several  scene  types,  including  a  downward-looking 
stare  scene  through  a  clear  atmosphere,  an  elevation  scan  over  varied 
background  including  a  feature  which  may  be  thin  high-altitude  cirrus,  and 
an  azimuthal  scan  over  a  background  scene  which  includes  broken  fair  weather 
cumulus  clouds  at  lower  altitudes. 

A. 1  Mountain  Terrain  Background  (Clear  Air) 

The  clear-air  comparisons  of  calculated  and  measured  radiance  levels 

'">8) 

are  made  for  a  staring  scene  identified  in  the  EAMM  Data  Report''  as 

Scene  Ic  and  Scene  Id.  This  scene  is  a  downward- looking  stare  (observer 

nadir  angle  *  1.8  -  3.1°)  at  the  Sacramento  Mountains.  The  background 

consists  of  forested,  hilly  terrain  with  the  elevation  varying  between  7000' 

and  8500'  (above  sea  level),  as  determined  from  U.S.  Geological  Survey 

topographic  maps.  For  the  three-minute  stare  sequences,  the  radiometer 

data  show  very  little  temporal  or  spatial  variation,  and  the  typical  scene 

2 

radiance  values  are  0.6  pW/cm  /sr/band  for  Filter  1  (wideband  filter)  and 
2  ( ic;) 

0.1  pW/cm  /sr/band  for  Filter  2  (narrow  band  filter).'' 

The  Aerodyne  background  radiance/cloud  model  was  used  to  calculate 
expected  levels  of  background  radiance  at  the  balloon  altitude  for  Scene  Ic 
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Murphy,  R.E.,  Cook,  F.H. ,  Grossbard,  N.G. ,  Grieder,  K.F.,  and  Yap,  B.K., 
"Balloon  Altitude  Mosaic  Measurements  (BAMM)  Data  Report  II,"  Air  Force 
Geophysics  Laboratory,  Hanscom  AFB,  MA  01731  (1979). 


and  Scene  Id.  A  baseline  calculation  using  the  "Holloman"  radiosonde  data 
(28  September  1978,  0630  MST)  was  performed.  Further  calculations  vtu-  au.de 
to  test  the  sensitivity  of  the  calculations  to  atmospheric  temperature  and 
humidity  profiles.  The  in-band  radiances  for  BAMM  Filters  1  and  2  were 
determined  for  each  case  by  integrating  the  calculated  spectral  radiance 
over  the  filter  bandpasses,  including  the  transmission  functions  of  the 
filters.  The  input  parameters  for  these  calculations  are  given  in  Table  2. 

A  uniform  surface  emissivity  of  0.9  and  reflectivity  of  0.1  were  assumed 
for  the  calculations.  A  high  emissivity  for  botanical  materials  is  expected 
due  to  their  high  cellulose  and  water  contents 


Table  2.  Baseline  Input  Parameters  for  Background  Radiance 
Calculations  (Scene  Ic  and  Scene  Id) 


Observer  altitude 

98,000  ft.  (29.9  km) 

Background  altitude 

7,000  ft.  (2.13  km) 

Observer  Angle  (zenith) 

2.8° 

Solar  angle  (zenith) 

48.3° 

Background  emissivity 

0.9  1  (assumed  constant 

Background  reflectivity 

0.1  )  with  wavelength) 

Background  temperature 

13°C 

(Holloman  radiosonde  data) 

- __i 

The  results  of  the  baseline  calculation  are  given  in  Table  3.  The 
calculated  in-band  radiances  are  shown  in  the  table,  along  with  measured 
in-band  radiances  from  Scene  I.  The  contributions  of  the  various  radiance 
components  (ground  thermal  emission,  ground  solar  scatter,  and  atmospheric 
thermal  emission)  to  the  total  calculated  values  are  also  given  in  the  table. 
__ 

Woli ,  W.L.,  and  Zissis,  G. J. ,  Editors,  The  Infrared  Handbook,  (Office  of 
Naval  Research,  Department  of  the  Navy,  Washington,  DC,  1978). 
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Table  3.  Baseline  Radiance  Calculation  for  Scene  1 


2 

In-Band  Radiance,  W/cm  /sr 

Filter 

Data 

Calculations 

Total 

Thermal 

Solar 

Scatter 

Atmospheric 

1 

2 

6  x  10“ 7 

1  x  10“7 

3.86  x  10“7 
0.70  x  10"7 

5.18  x  10“8 
0.76  x  10“8 

2.42  x  10"7 

0.38  x  10“7 

9.22  x  10“8 

i 

2.44  x  10~8 
_ 

The  sensitivity  of  the  calculations  to  changes  in  the  atmospheric  temper¬ 
ature  and  humidity  profiles  is  indicated  by  the  results  given  in  Table  4.  The 
calculated  in-band  radiances  are  shown  in  this  table  for  various  atmospheric 
profiles,  otherwise  using  the  baseline  input  of  Table  2.  Looking  at  the  ratio 
of  calculated  Filter  1  and  Filter  2  band  radiances  in  Table  4,  we  see  that  the 
ratio  of  the  two  bands  is  fairly  insensitive  to  the  atmospheric  profile. 
Reasonably  close  agreement  is  obtained  among  the  calculations  made  with  the 
two  sets  of  radiosonde  data  and  the  midlatitude  summer  model.  The  "Jallen" 
radiosonde  data  correspond  to  an  ascension  of  time  of  1330  MST  on  Sept.  28, 
whereas  the  "Holloman"  ascension  time  was  0630  MST.  The  atmospheric  profile 
for  the  Jallen  (early  afternoon)  data  is  warmer  and  drier,  resulting  in 
reduced  atmospheric  absorption.  This  effect  causes  significant  increases  in 
reflected  solar  energy  as  measured  at  the  observer  altitude.  The  proportion 
of  the  total  radiance  which  is  due  to  reflected  sunshine  also  increases  due  to 
the  double  transit  of  this  radiation  through  the  atmosphere.  The  fractions 
of  the  total  radiance  which  are  due  to  reflected  solar  radiation  are  given  in 
Table  5  for  the  two  atmospheric  profiles. 


Table  4.  Calculated  In-Band  Radiance  for  Several 
Atmospheric  Profiles 


Atmospheric  Model 

r-  — 

Radiance, 

W/cm^/sr 

Nl/N2 

Filter  1 

Filter  2 

Holloman  Radiosonde 

0.39  x  10"6 

0.70  x  10-7 

5.57 

Jallen  Radiosonde 

0.52  x  10‘6 

0.99  x  10-7 

5.26 

Midlatitude  Summer 

0.48  x  10-6 

0.84  x  10-7 

5.74 

1962  U.S.  Standard 

0.89  x  10-6 

r-~ 

1 

o 

rH 

X 

oo 

• 

rH 

4.83 

Observed 

0.6  x  10  6 

1.0  x  10"7 

6.0 

Table  5.  Fraction  of  Calculated  In-Band  Radiance  Due 
to  Reflected  Solar  Radiation 


1 

Reflected  Solar  Energy  j 

Atmospheric  Model 

(Fraction  of  Total  Radiance)  ; 

Filter  1 

Filter  2 

Holloman  Radiosonde 

0.54 

Jallen  Radiosonde 

0.64 

These  calculations  indicate  that  the  modified  LOWTRAN  atmospheric 
transmission  and  background  radiance  model  predicts  levels  of  background 
radiance  for  this  scene  which  are  in  good  agreement  with  the  BAMM  radiometer 
measurements.  The  predicted  signal  ratios  between  the  two  filters  are 
consistent  with  the  data.  The  absolute  predicted  radiance  levels  can  be 
shifted  somewhat,  depending  on  the  choice  of  background  emissive  and  reflective 
properties.  In  these  calculations,  approximately  50-70%  of  the  total  back- 
bround  radiance  is  due  to  reflected  sunlight.  This  fraction  is  highly 
dependent  on  the  background  elevation  assumed  in  the  calculations.  For  lower 
background  elevations  the  atmospheric  transmission  decreases  significantly. 
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primarily  because  of  enhanced  absorption  by  water  vapor,  and  the  fractional 
contribution  of  reflected  sunlight  to  the  total  background  radiance  would  be 
smaller.  The  model  calculations  discussed  in  the  next  subsection  illustrate 
this  effect. 


4.2  Elevation  Scan:  Desert /Mountain  Background  With  Possible  Cirrus  Cloud 

The  first  portion  of  the  Scene  I  data  (Scene  la  and  Scene  lb)  discussed 
(28) 

in  the  BAMM  report  consists  of  an  elevation  scan,  beginning  at  a  maximum 
nadir  angle  of  approximately  69°  and  scanning  down  to  a  minimum  angle  of 
approximately  1.8°  (downward  viewing).  A  radiometer  (wideband  filter) 
radiance-time  plot  for  this  scene  is  shown  in  Fig.  12. 

The  plots  of  absolute  radiance  versus  time  for  the  Scene  I  elevation 
scan  (Scene  la  and  Scene  lb)  show  several  characteristic  features.  A  general 
feature  is  the  increase  in  average  signal  level  as  the  elevation  angle 
decreases.  The  mean  signal  level  at  maximum  nadir  angle  (69.47°)  is 
approximately  55%  of  its  value  at  the  minimum  nadir  angle.  The  increase  in 
signal  with  decreasing  nadir  angle  is  mainly  due  to  the  increase  in  atmos¬ 
pheric  transmittance  for  a  shorter  path  from  the  ground  to  the  balloon. 
Calculated  atmospheric  spectral  transmittance  for  observer  nadir  angles  of  0° 
and  60°  are  shown  in  Fig.  13.  (A  background  elevation  of  4100  ft  is  used  in 
the  calculations.)  These  LOWTRAN  calculations,  which  include  only  atmospheric 
and  surface  thermal  emission,  show  that  the  expected  signal  level  at  a 
nadir  angle  of  69°  is  approximately  75%  of  its  value  for  downward-looking. 

The  measured  signal  also  includes  a  solar  scattering  component.  This 
component  becomes  more  strongly  attenuated  at  the  high  nadir  angles  because 
the  solar  radiation  which  is  reflected  makes  a  double  transit  through  the 
atmosphere.  If  diffuse  solar  scattering  with  a  terrain  reflectance  of  0.1 
is  assumed,  the  LOWTRAN  model  predicts  a  variation  in  absolute  radiance  level 
which  closely  matches  the  measured  dependence  on  nadir  angle.  This  prediction 
is  the  dotted  curve  labeled  "Aerodyne  calculations"  in  Fig.  12.  The  fractional 
contribution  of  the  solar  scatter  component  increases  from  approximately  20° 


Comparison  of  Background  Radiance  Calculations  With  BAMM 
Radiometer  Measurements  for  Scene  la  and  Ib(24) 
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Figure  13. 
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ac  70°  to  37%  at  0°.  This  fractional  contribution  is  less  than  the  scattering 
contributions  discussed  in  the  previous  comparison  (stare  scene)  because  the 
earth  background  was  at  a  lower  elevation. 

In  addition  to  the  overall  increase  in  signal  level  which  is  due  to 
the  Increasing  atmospheric  transmittance,  two  other  features  of  the  radiance¬ 
time  plots  for  Scene  la  and  Scene  lb  are  noteworthy.  The  first  of  these  is 
the  large  increase  in  signal  which  occurs  in  Scene  la  at  approximately 
37330  sec  (10:22:10).  The  signal  in  all  channels  increases  by  nearly  a 
factor  of  two.  This  signal  increase  is  correlated  with  a  blurring  of  the 
scene  in  the  television  pictures.  This  blur,  which  does  not  totally 
obscure  the  terrain  background,  has  not  been  positively  identified.  A  likely 
explanation  is  that  it  is  a  thin  cirrus  cloud,  however.  This  possibility  is 
examined  in  more  detail  below.  In  the  latter  part  of  the  elevation  scan 
(Scene  lb),  the  signal  in  all  channels  shows  variations  which  are  correlated 
with  spatial  variations  in  the  background.  The  background  scene  in  this 
portion  of  the  scan  is  the  Sacramento  Mountains  with  sunlit  peaks  and 
shadowed  regions  where  the  solar  scatter  contribution  to  the  signal  is 
reduced.  The  measured  signal  level  is  modulated  by  scanning  over  this 
spatially  structured  terrain. 

The  peak  in  the  infrared  signal  in  Scene  la,  which  is  correlated  with  a 
blur  in  the  television  pictures,  suggests  the  possibility  of  a  thin  high 
altitude  cloud.  Such  a  cloud  can  transmit  much  of  the  upwelling  radiance 
from  the  earth  surface  and  the  lower  part  of  the  atmosphere  in  addition  to 
reflecting  solar  radiation.  In  particular,  radiation  in  the  2.7  vim  water  vapor 
bands  which  is  scattered  by  high  altitude  clouds  can  lead  to  an  enhanced  back¬ 
ground  signal  because  the  cloud  is  above  most  of  the  atmospheric  water 
vapor  attenuation.  The  two-stream  cloud  radiance  model  developed  in  Section  2 
has  been  used  to  illustrate  these  effects  and  to  examine  the  possibility  that 
the  signal  deriving  from  the  "blur"  could  be  caused  by  a  cirrus  cloud  at  high 
altitudes.  The  results  of  these  calculations  are  summarized  in  Table  6, 
where  the  calculated  in-band  radiance  (BAMM  Filter  1)  for  a  background 
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Table  6.  Calculated  In-Band  Radiance  (BAMM  Filter  1)  For 
a  Background  Scene  Including  An  Ice  Cloud  at 
50,000  ft  Altitude. 


CLOUD  THICKNESS 

(kn) 

TRANSMITTANCE 

(at  *3.2  ^im) 

SCENE  RADIANCE 
CfxV/  ctn2/sr) 

0 

1.0 

0.6 

0.1 

0.97 

1.1 

0.2 

0.95 

1.4 

1.0 

0.82 

3.6 

Opaque  Cloud 

0. 

14. 

i 

including  a  cloud  at  50,000  ft  altitude  is  shown  for  a  range  of  cloud  thick- 

_3 

nesses.  A  particle  concentration  of  0.01  cm  was  assumed  for  the  calculations 
Scattering  and  absorption  properties  for  a  distribution  of  ice  particles  with 
a  mean  particle  radius  of  32  pm  were  used  in  the  calculations. 

The  results  in  Table  6  demonstrate  that  scattering  of  sunlight  from 
thin  high-altitude  clouds  can  result  in  significantly  enhanced  background 
signals,  when  the  scattering  effects  are  not  blurred  by  atmospheric  attenua¬ 
tion.  The  calculated  band  radiance  for  BAMM  Filter  1  for  a  cloud  with  a 

transmittance  of  0.95  at  3.2  pm  (0.2  km  thickness)  is  more  than  twice  the 

-3 

clear-air  signal.  The  particle  concentration  of  0.01  cm  is  consistent 
with  or  even  smaller  than  typical  particle  concentrations  in  cirrus  clouds. 

The  predicted  signal  enhancement  for  such  a  cloud  is  consistent  with  the 
observed  enhancement  evident  in  the  radiometer  data  of  Fig.  12. 


The  cloud  altitude  has  a  strong  effect  on  possible  background  signal 
enhancement  due  to  solar  scattering.  The  altitude  effect  is  illustrated  in 
Fig.  14,  which  shows  a  comparison  of  calculated  cloud  radiance  with  the 
clear-air  background  radiance  for  two  cloud  altitudes.  At  high  altitudes 


(12.2  km),  the  cloud  Is  above  most  of  the  atmospheric  water  vapor 
attenuation  and  the  spectral  radiance  In  the  2. 5-2. 8  pm  region  is  greatly 
enhanced  by  the  presence  of  the  cloud.  Atmospheric  attenuation  by  CC>2  is 
evident  in  the  calculated  cloud  spectral  radiance.  At  lower  altitudes 
(4.57  km),  the  cloud  signal  is  strongly  attenuated  below  2.85  pm 
by  atmospheric  H^O.  Beyond  2.85  pm,  the  cloud  signal  may  be  either  larger  or 
smaller  than  the  clear-air  background  signal,  depending  on  the  efficiency  of 
the  directional  scattering  of  the  cloud  compared  with  diffuse  scattering 
from  the  terrain  background. 

The  BAMM  television  pictures  offer  no  indication  of  the  altitude  of  the 
"blur"  which  is  evident  in  Scene  la.  As  indicated  in  the  discussion  above, 
the  infrared  data  for  the  scene  are  consistent  with  model  predictions  for 
an  optically  thin  cloud  at  high  altitudes.  The  cloud  model  was  also  used 
to  simulate  ground  fog  for  the  conditions  of  Scene  la,  by  locating  the 
model  "cloud"  as  a  layer  just  above  ground  elevation.  For  this  case,  the 
model  calculations  show  that  the  fog  inhibits  transmission  at  wavelengths 
where  molecular  attenuation  is  not  important,  and  the  layer  scatters  sun¬ 
light  less  efficiently  than  the  terrain  background.  Strong  H^O  and  C02 
attenuation  prevents  any  significant  contribution  to  the  total  signal  below 
2.85  pm.  The  net  effect  is  that  the  predicted  radiance  from  the  background 
with  ground  fog  is  reduced  compared  to  the  clear-air  terrain  background, 
toith  a  maximum  signal  reduction  of  25%  for  optically  thick  fog. 

The  model  predictions  discussed  in  this  section  show  that  high  altitude 
clouds  can  produce  strong  IR  background  radiance  signals  in  the  2. 7-3. 2  urn 
region,  and  that  the  signals  may  be  significantly  larger  than  the  background 
radiance  for  clear-air  conditions.  The  effect  is  significant  even  for  an 
optically  thin  cloud.  This  effect  has  probably  been  observed  in  BAMM  Scene  la; 
however  the  television  pictures  do  not  provide  conclusive  proof. 
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4.3  Cumulus  Clouds  (Low  Altitude) 

The  data  from  the  BAMM  Holloman  radiometer  test  flight  show  significant 

background  radiance  variations  in  Scene  II,  which  is  an  azimuthal  scan  over 

(28) 

a  mountain  background,  with  broken  cumulus  clouds.  '  The  significant  back¬ 
ground  features  are  sunlit  terrain,  shadowed  terrain  (from  cloud  shadowing), 
and  cumulus  clouds  which  are  optically  thick  in  the  visible.  The  individual 
clouds  are  typically  of  a  size  such  that  they  fill  most  of  the  field  of  view 
of  one  radiometer  pixel.  A  radiance  time  profile  for  BAMM  Scene  lib  (Fig.  29 
reproduced  from  Ref.  (28))  is  shown  in  Fig.  15.  The  plot  shows  the  measured 
radiance  in  one  pixel  for  the  three-minute  scan,  using  Filter  1  (the  wideband 

filter).  By  comparison  of  the  data  for  all  of  the  detector  elements,  the 

2 

maximum  measured  radiance  for  the  scan  is  approximately  1.0  pW/cm  /sr/band, 

2 

with  minimum  radiance  levels  of  approximately  0.45  pW/cm  /sr/band.  A  similar 

scan  in  Scene  11a  with  Filter  2  (narrow  filter)  resulted  in  measured  radiance 

2 

variations  from  approximately  0.088  pW/cm  /sr/band  to  maximum  values  of 

2 

approximately  0.16  pW/cm  /sr/band.  Comparison  of  the  radiance  time  profiles 
with  the  television  pictures  indicates  that  the  largest  signals  correspond  to 
terrain  backgrounds,  intermediate  signal  levels  correspond  to  clouds  or  mixed 
cloud-terrain  backgrounds,  and  the  lowest  signal  levels  correspond  to  cloud 
shadows . 


The  two-stream  cloud  scattering  model  was  used  to  calculate  the  expected 
radiance  levels  for  clouds,  cloud  shadows,  and  sunlit  terrain  for  the  conditions 
of  the  Scene  II  azimuth  scans.  The  background  elevation  was  determined  from 
USGS  topographic  maps,  and  cloud  heights  were  estimated  from  the  television 
pictures  (based  on  the  relative  cloud  and  cloud  shadow  positions).  The  cloud 
calculations  were  carried  out  for  the  range  of  relative  observer  and  solar 
azimuth  angles  corresponding  to  the  measurement  scan,  however  the  results 
were  Insensitive  to  azimuth  angle.  The  input  conditions  for  the  calculations 
are  summarized  in  Table  7. 


Table  7.  Input  Parameters  for  Low  Altitude  Cumulus 
Cloud  Calculations  (BAMM  Scene  II) 


Terrain  Background: 

Elevation 

Temperature 

Emissivity 

Reflectivity 

Atmospheric  Profile 

Cloud  Background: 

Cloud  Top  Elevation 
Cloud  Thickness 
Particle  Density  (Liquid  H^O) 
Mean  Particle  Radius 

Observer  Geometry: 

Elevation 
Azimuth  Angle 
Nadir  Angle 

Solar  Angles: 

Zenith 

Azimuth 


2.59  km 
20°C 
0.90 
0.10 

"Holloman 


5.2  km 
1  km 

500  cm-3 
4  pm 

30.2  km 
148.4-247 
24.4° 

38.4° 

153° 


constant  with 
wavelength 

'  radiosonde 
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The  calculated  in-band  radiance  levels  for  the  various  background 
features  are  shown  in  Table  8.  The  expected  radiance  levels  for  cloud 
shadows  are  calculated  by  setting  the  solar  scatter  contribution  to  zero.  T 
calculation  thus  includes  only  thermal  radiation  from  the  terrain  background 
and  from  the  atmosphere,  with  atmospheric  attenuation.  The  predicted  cloud 
shadow  signal  levels  for  both  filters  are  in  excellent  agreement  with  the 
measured  minimum  signals  during  the  scans. 


Table  8.  Calculated  Radiance  Levels  (In-Band)  for  Background 
Features  in  BAMM  Azimuth  Scans  (Scene  II) 


Background 

2 

Radiance,  pW/cm  /sr/band 

Filter  1  (Wide  Filter) 

Filter  2  (Narrow  Filter) 

Sunlit  Terrain 

2.0 

0.39 

Cumulus  Clouds 

0.64 

0.13 

Cloud  Shadows 

0.41 

0.091 

An  apparent  discrepancy  between  the  measurements  and  the  calculations 
arise  in  comparing  the  predicted  radiance  levels  for  sunlit  terrain  with  the 
maximum  measured  signal.  The  predicted  sunlit  terrain  radiance  is  a  factor 
of  2  larger  than  the  highest  observed  signals.  In  the  case  of  terrain  back¬ 
ground,  the  solar  scatter  term  is  calculated  by  using  a  specified  input 
reflectivity  of  0.10.  Variations  of  this  quantity  with  wavelength,  as  well 
as  a  lower  effective  reflectivity  due  to  surface  roughness  (terrain  self¬ 
shadowing)  may  account  for  overprediction  of  the  solar  scatter.  In  addition, 
the  background  elevation  assumed  for  the  calculation;-  (2.59  km  or  8500')  r.iv 
be  higher  than  the  mean  elevation  during  the  scan,  in  which  case  the  atr-esphe 
attenuation  is  underestimated. 


The  qualitative  trend  of  highest  IR  signals  being  correlated  with  sunlit 

terrain,  intermediate  signal  levels  being  correlated  with  clouds,  and  minimum 
signal  levels  being  due  to  cloud  shadows  is  predicted  by  the  cloud  model.  Ex 


correlation  of  the  television  pictures  with  the  infrared  data  is  difficult, 
however,  because  the  instantaneous  field  of  view  for  any  detector  element 
typically  Includes  some  combination  of  cloud,  cloud  shadow,  and  sunlit 
background.  The  experimental  differences  between  cloud  radiance  and  sunlit 
terrain  radiance  for  this  scene  are  therefore  difficult  to  assess. 
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Scene  III  from  the  BAMM  radiometer  test  flight v  ’  is  a  staring  scene 

at  a  low  altitude  cumulus  cloud.  During  the  course  of  the  stare  sequence, 

the  combination  of  cloud  and  platform  drift  causes  the  f ield-of-view  for 

several  of  the  detector  elements  to  shift  from  mostly  cloudy  to  mostly  clear. 

Radiance-time  profiles  for  this  sequence  show  negligible  long-term  drift  of 

the  in-band  radiance  signal,  suggesting  little  difference  between  cloud  and 

terrain  backgrounds  in  this  instance.  Figure  16  shows  the  measured  radiance- 
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time  profil  for  Detector  7.  By  comparison  with  the  television  pictures, 
the  field-of-view  for  this  detector  changes  from  completely  cloudy  (visibly 
thick  cloud)  to  almost  completely  clear  (in  the  visible)  over  the  course  of 
the  scan. 

The  experimental  results  from  the  radiometer  test  flight  suggest  that 
thick  cumulus  clouds  at  elevations  around  4.6  km  (15,000  ft)  have  little 
effect  on  the  band  radiance  in  the  2.7-3. 1  urn  region.  The  cloud  model 
calculations  predict  a  more  significant  effect,  however  the  strong  radiance 
variations  that  were  evident  for  the  high  altitude  cirrus  clouds  are 
moderated  by  atmospheric  attenuation  at  the  lower  altitudes.  Further 
comparisons,  especially  emphasizing  the  spectral  radiance  from  interferometer 
data,  are  required  to  validate  the  cloud  model  for  low  altitude  clouds. 
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5.  SUMMARY 


A  model  which  describes  the  scattering,  absorption,  and  emission  of 
radiation  by  atmospheric  water  or  ice  clouds  has  been  developed,  based  on  an 
approximate  solution  to  the  radiative  transfer  equation.  The  cloud  model  is 
combined  with  the  LOWTRAN  atmospheric  transmission  model  and  can  be  used  to 
predict  the  apparent  spectral  radiance  from  clouds  as  detected  by  a  downward- 
looking  observer,  including  the  effects  of  solar  scatter,  thermal  emission, 
transmitted  upwelling  radiation,  and  atmospheric  attenuation.  The  model  is 
completely  general  in  terms  of  sun-cloud-observer  geometry  and  applies  for 
clouds  of  arbitrary  optical  thickness. 

Calculations  using  the  model  have  been  compared  with  infrared  cloud 
data  from  several  sources.  Comparisons  with  spectral  data  in  the  2. 7-2. 9  pm 
and  2. 5-6.0  pm  regions  show  generally  good  agreement  in  the  spectral  dependence 
of  the  cloud  signal,  although  discrepancies  on  the  order  of  a  factor  of  two 
are  indicated  in  the  directional  dependence  (i.e.,  for  solar  scattering  angles 
of  10°  or  less).  More  detailed  comparisons  with  a  larger  set  of  data  are  required 
td  assess  the  practical  importance  of  the  various  approximations  in  the  model. 

Comparisons  of  model  predictions  with  data  from  the  BAMM  radiometer  test- 
flight  show  generally  good  agreement  in  terms  of  predicted  and  measured 
terrain  background  signals.  Both  the  data  and  the  calculations  suggest  that 
low  altitude  cumulus  clouds  are  not  an  important  source  of  background  signal 
variations  in  the  2. 7-3.1  pm  region.  More  detailed  data  comparisons,  parti¬ 
cularly  emphasizing  spectral  radiance,  are  required  to  establish  the  relative 
magnitudes  of  low  altitude  cloud  signals  compared  with  the  background  terrain. 

The  model  calculations  illustrate  that  high  altitude  cirrus  clouds  can 
cause  large  variations  in  the  upwelling  radiation  signal  in  the  2. 7-3.1  pm 


region.  The  band  radiance  for  background  scenes  including  such  a  cloud 
can  be  several  times  the  clear  air  radiance,  even  for  an  optically  thin 
cloud.  This  effect  arises  because  the  cloud  is  above  most  of  the  atmospheric 
water  vapor,  and  attenuation  of  the  scattered  sunlight  is  not  as  severe  as  in 
the  case  of  low  altitude  clouds.  An  analysis  of  elevation  scan  data  show  that 
solar  scattering  from  a  high  altitude  cirrus  cloud  is  a  plausible  explanation 
for  the  large  IR  signal  increase  associated  with  apparent  haze  in  HAMM  Scene 
la.  A  consideration  of  potential  signal  variations  from  such  clouds  is  likely 
to  be  important  in  the  design  of  downward-looking  surveillance  systems. 

The  cloud  model  discussed  in  this  report  successfully  describes  the 
background  IR  signals  from  different  scenes.  Including  both  thin  and  thick 
clouds.  The  model  provides  a  sueful  tool  for  analyzing  and  interpretating 
background  scenes  observed  from  the  BAMM  platform.  These  analyses  are 
important  in  assessing  the  upwelling  radiances  which  would  be  measured 
from  a  satellite  platform. 
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